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Abstract. Clebsch-Gordan coefficients of SU(2) and SU(1,1) are defined as 
eigenfunctions of a linear operator acting on the tensor product of the Hilbert spaces for 
two irreps of these groups. The shifted harmonic approximation is then used to solve 
these equations in asymptotic limits in which these eigenfunctions approach harmonic 
' oscillator wave functions and thereby derive asymptotic expressions for these Clebsch- 

Gordan coefficients. 
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1. Introduction 



The shifted harmonic approximation (SHA) is an approximation to a technique for 
realizing a set of operators, defined initially as linear transformations of a finite vector 
space, as differential operators. It was introduced for the purpose of understanding the 
nature of phase transitions in systems with su(2) or su(2) © su(2) spectrum generating 
algebras [1] and subsequently applied to a model with an su(l,l) © su(l,l) spectrum 
generating algebra [2]. More recently it has been developed for application to a multi- 
level pairing model [3], for which the spectrum generating algebra is a direct sum of 
multiple su(2) algebras. It can be regarded as a technique for contracting a Lie algebra 
representation to that of a simpler algebra in well-defined limiting situations. 

In this paper, we show how the SHA provides approximate expressions for SU(2) 
and SU(1,1) Clebsch-Gordan coefficients that become precise in certain asymptotic 
limits in which they approach harmonic oscillator wave functions when regarded as 
functions of appropriate parameters. Moreover, we show by examples, that these limits 
are approached very rapidly and provide remarkably accurate approximations more 
generally. 

In addition to the examples given above, the groups SU(2) and SU(1,1) and the 
coupling of their irreducible representations have numerous applications in quantum 
optics In this context, asymptotic SU(2) and SU(1,1) CG coefficients are of huge 

| A highly non-exhaustive list contains, for instance [5]. 
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potential value in problems in which a prohibitive amount of time is often needed to 
compute the required coefficients numerically from the known analytical formulae. In 
contrast, the asymptotic coefficients given here can be calculated in fractions of a second. 
Asymptotic SU(2) CG coefficients have also been connected to a tight-binding model of 
a one-dimensional potential [6]. 

Clesbch-Gordan (CG) coefficients (si,M — m, s 2 ,m\SM) for the group SU(2) are 
required when two systems of spin si and s 2 are coupled to total spin S: 

\SM) = ^2\s 2 m) <g> \si,M -m)(s 1 ,M -m, s 2 ,m\SM). (1) 

m 

In Eqn. (JIJ), the projections along the common quantization axis of the spins s±, s 2 , and 
S are M — m,m and M, respectively. Eqn. is also applicable to a wide variety of 
other quasi-spin systems having states that carry su(2) irreps. 

Because of the connection between SU(2) and SO (3), CG coefficients also appear 
in problems where products of spherical harmonics and other related special functions 
occur naturally. Furthermore, the ubiquity of su(2) as a subalgebra of other Lie algebras 
makes the SU(2) CG coefficients (and those of SU(1,1) for non-compact groups) an 
important ingredient in the computation of coupling coefficients for higher groups. 

Much is known about SU(2) CG coefficients: any coefficient can be obtained in 
closed form using an expression containing a sum of square roots of rational factors. 
However, for small values of m and M, the complexity of this sum increases rapidly with 
S, and numerous asymptotic estimates for large S have been developed to understand 
and quickly evaluate CG coefficients in this regime. 

A first consideration of asymptotic limits for SU(2) coefficients is due to Wigner 
[7J . He approached the problem from a semi-classical perspective and obtained average 
expressions for the coefficients when all momenta were large; unfortunately, Wigner's 
result did not capture the essentially oscillatory nature of asymptotic CG coefficients 
and left the door open for further studies. Brussaard and Toelhoek [8] used the WKB 
method to connect reduced SU(2)-Wigner functions with asymptotic CG coefficients. 
Their result, which applies to an asymptotic limit in which just two angular momenta Si 
and S are large, significantly improved on Wigner's because it correctly gave the sign of 
the coefficient and turns out to be reasonably accurate, even for modest values of Si and 
S. Finally, Ponzano and Regge [9], in their authoritative work grounded almost entirely 
on geometrical arguments, expanded on [7] and [8] by providing accurate expressions 
for CG coefficients valid for three large angular momenta. 

The insights provided by [Hj and [9] have been a source of inspiration for many 
subsequent authors. A detailed survey of the literature pre-1988 can be found in [TO] . 
Comparatively recent work, limited to asymptotic SU(2) CG coefficients with all three 
momenta large, include, for instance, Refs. [HI [12], wherein a detailed and systematic 
review of the results of Ponzano and Regge can be found, with emphasis on closing loose 
ends in some of their "heuristic" arguments. 

In this paper we show that for asymptotically large values of s% and s 2 , and for 
finite values of n and M, the SU(2) CG coefficients rapidly approach the asymptotic 
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expressions 

(si,M— m, s 2 m\si + s 2 — n, M) 



-y H n (a(m - x ))e-^ m -^\ (2) 



. v / 7r2 n n! , 

where H n is a Hermite polynomial, with parameters defined by a simple SHA algorithm. 
It is also shown that these asymptotic results have analytic approximations which remain 
precise in the asymptotic limit but are approached somewhat less rapidly; they are given 
by the explicitly expression 

(simi s 2 m 2 |si + S2- n, M) 

~ (-!)•( °)* g .(„ (3) 
V yj 7i2 nl / V o~i + cr 2 ' 

with 

a 4 = + = 

o~i = a/si(si + 1) , a 2 = y/s 2 (s 2 + 1) • (5) 

Basis states for unitary irreps of the positive discrete series of su(l,l) within the 
tensor product of two such irreps are also given by linear combinations 

\Kn) = \hn 2 } <8> IhnJ (hk 2 Nm\Kn), (6) 

nin2 

where (kik 2 Nm\Kn) is an SU(1,1) CG coefficient related to a non-vanishing coefficient 
in the notation of Van der Jeugt [13], by 

(kik 2 Nm\Kn) — C(k\, ni, k 2 , n 2 ; K, n) (7) 

with N = ni + n 2 and m = n 2 — ri\. Note that an SU(1,1) CG coefficient vanishes unless 
K + n = k\ + n\ + k 2 + n 2 = k\ + k 2 + N. The notation is clarified in Sectj3j 

Various closed form expressions for the SU(1,1) coefficients can be found in the 
literature. The exact expressions of immediate relevance to our work are given, for 
instance, in [13JQ3]. In addition, several other authors [15] have published expressions 
spanning a variety of couplings of representations (not only of two positive discrete 
series) and a number of different bases. 

In spite of these exact results, there is limited knowledge of the asymptotic behavior 
of the coefficients, although some basic results valid for the coupling of two irreps of the 
positive discrete series with k\ — > oo and k 2 finite can be found in [16], a result that 
proves to be useful in constructing states in odd nuclei p2] within the context of the 
nuclear symplectic model. 

The SHA approach is applied in Sect. [3] to derive asymptotic expressions for the 
SU(1,1) CG coefficients of Eqn.flTj) in an altogether different regime to that studied in 
[T7] . The result, similar to that given above for a class of SU(2) coefficients, is that 

(k x k 2 Nm\Kn) 

- (-l)^ +m (^^)^n(a(m-x ))e^ a2 ^ 2 . (8) 
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when N = n 1 +n 2 becomes asymptotical large and n remains finite, and with parameters 
defined by a simple SHA algorithm given in Sect. 13.11 It is also shown that precise 
asymptotic expressions which, however, are approached somewhat less rapidly, are given 
by Eqn. (J8]) with the explicit parameter values 

2N^T 2 ' X ° ~ N + 4 V /^' [ ' 

where 

N , N 

«i = «i + — , k 2 = k 2 + — . (10) 

The SHA method described in this paper complements other approaches where 
the asymptotic behavior is examined using WKB methods. As mentioned briefly in 
the Discussion, the SHA technique can be regarded as a procedure for obtaining the 
contraction of a Lie algebra to a simpler Lie algebra that is appropriate in certain 
limiting situations. Such contractions are known to lead to valuable insights and useful 
approximation procedures in physics as evidenced, for example, in the approach of 
quantum mechanics to classical mechanics as h or some other scale parameter approaches 
zero, in the approximation of fermion pair algebras by boson algebras in the random- 
phase approximation of many-body theory, and in the bosonic behaviour of large atomic 
samples [T8] . 



2. Asymptotic SU(2) Clebsch-Gordan coefficients 

Let {S + , S_, S } satisfy the usual commutation relations of the complex extension of 
the su(2) Lie algebra: 

[S ,S±] = ±S ± , [S+,S-]=2S . (11) 

Basis states, \sm), for an su(2) irrep are then defined by the equations 

So\sm) = m\sm), (12) 

S±\sm) = y/ (s =F rn)(s ±m + 1) \s, m ± 1). (13) 

Coupled basis states for irreps of su(2) within the tensor product of two such irreps 
are given by 

\SM) = \ s 2 m 2) ® \s 1 m 1 )(s 1 m 1 s 2 m 2 \SM), (14) 

m\m2 

where (s\m\ s 2 m 2 \SM) is an SU(2) CG coefficient. These coefficients are equal to the 
overlaps of coupled and uncoupled tensor product states 

(simi s 2 m 2 \SM) = [(simi| <8> (s 2 m 2 \}\SM). (15) 

For convenience, we denote the tensor product states by 

\s 1 s 2 Mm) = \s 2 m 2 ) ® (16) 

with M = mi + m2 and m = m 2 . The SU(2) CG coefficients are then the overlaps 

(sxmi s 2 m 2 \SM) = (s 1 s 2 Mm\SM) = [(sm^ g) (s 2 m 2 \}\SM). (17) 
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They are defined (to within arbitrary phase factors) by the requirement that the states 
{I^M)} satisfy the eigenvalue equations 

S Q \SM) = M\SM), (18) 
S + S.\SM) = [S(S + 1) - M[M - 1)] \SM), (19) 

with 

S = Si + Si S ± = Si + S 2 ± . (20) 
2.1. The shifted harmonic approximation 

We now determine the CG coefficients using the SHA and show them to be precise in 
the limit of asymptotically large values of si, s 2 , and S and finite values of M and 
n = si + s 2 - S. 

The desired CG coefficients are first expressed as overlap functions, in the form 

i> s n 1S2M (m) = (s 1 s 2 Mm\s l + S2-n,M). (21) 

This notation is introduced with the intention that, for given values of si, s 2 , M and n, a 
set of CG coefficients can be regarded as the values of a function, ip^ S2M , of the discrete 
variable m. Moreover, the state \SM) with S = si + S2—n is completely determined by 
the values, {ip^ lS2M (m),m = —s 2 , ■ ■ ■ , s 2 }, which implies that ip^ lS2M can be interpreted 
as a wave function for the state \si + s 2 — n, M). 

The operators So, S± are now mapped to operators on these wave functions defined 

by 

5„C 1S2M H = (sis 2 Mm\S u \ Sl + s 2 -n,M) 

= Y,( s i s 2Mm\S u \s 1 s 2 M'p}^ S2M '(p). (22) 

M'p 

Therefore, because the state \SM) is an eigenstate of an d S + S-, the function ^ lS2Af 
should likewise be an eigenfunction of So and S+S- with the same eigenvalues. With 
the expansions of Eqs. f|T6l) and fl20l . we obtain 

Sj_^M(m) = f (m)r n lS2M M 

+ fi{m)r n 1S2M (m + 1) + U{m)r n XS2M {m - 1), (23) 

where 

foijn) = (sis 2 Mm\S + S~\sis 2 Mm) , 
/i(m) = (sis 2 Mm|£+£_|sis 2 M,m + l), 

f-iim) = (s lS2 Mm|^ + ^_|siS 2 M,m- 1) = f^m - 1). (24) 
Equation ( 1231) can now be expressed in terms of finite difference operators, defined 

by 

Ai[)(m) = l(ip(m + 1) -ij)(m- I)), (25) 
A 2 ^(m) = ij(m+ 1) -2xjj(m) +ip{m- 1), (26) 
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with the result that 

A(m)^(m + 1) = Mm) (l + A + ±A 2 ) ^(m), (27) 

/i(m - lMm - 1) = (l - A + ±A 2 ) [AH^H], (28) 
Eqn. f )23|) then becomes 



5 +l S_C lS2M (m) = / (m)C 1S2M (m) + /i(m)[l + A + ±A 2 ]^ lS2M < 



m 



+ [1 - A + ±A 2 ] [/i(m)< 1S2M (m)] , (29) 
and gathering terms leads to the expression of S + S- as the difference operator 

= F(m) + AA(m)A, (30) 

where 

F(m) = / (m) + A(m) + /x(m - 1). (31) 
With = (Si + S 2 )(Si + S 2 ), we also determine that 

f Q (m) = a\ — (M — m)(M - m - 1) + a\ - m(m - 1), (32) 
/i(m) = ^[a 2 - (M — m)(M — m — 1)] [cr| — m(m + 1)], (33) 

where 

^ = 8^ + 1), 2 = 1,2. (34) 

To determine asymptotic expressions for the functions, ipn lS2N , as eigenf unctions of 
we now make the continuous variable approximation of extending these functions 
of the discrete variable m to functions, \l/^ lS27V , of a continuous variable x with the 
property that 

K 1S2N (x)=^ n 1S2N (x), (35) 

whenever a; is in the domain of the discrete variable m. In this approximation, which is 
valid in the asymptotic limits in which ty^ S2N becomes a smooth function, the difference 
operators can be replaced by differential operators: 

The functions F, /o and fi of Eqns. (13T]) - (133|) are similarly extended to the continuous 
variable x, and the operator S + S_ becomes the differential operator 

<S+<S„ -> F(x) + Dfi(x)D. (37) 

Provided the extension of ip^ lS2M (m) to the function \I/^ lS2M (:r) does not require the 
latter to be non-zero for any x that is outside of the limits for m, it is seen that fi(x) 
is real for all x for which ^! s ^ S2M (x) is non-zero. The limits on the values of m are seen, 
from Eqn. (TTB1) . to be such that — s 2 < m < s 2 and —si < M — m < S\. We will denote 
the upper and lower limits on the value of m, by m max and m min , respectively. Then, 
because the norm of the function ip^ lS2M is given by 

Hn 1S2M \\ 2 = E Wn S2M (rn)\ 2 , (38) 

771 — Tn-min 
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it follows that the corresponding smooth function \P* 1S2M (x) should have norm given by 

rrrimax 

U^ 2 M||2 = / \ss»^ M { x )\ 2 dx. (39) 



It is also seen that, when \E^ lS2M (x) is zero for all x > m max and all x < m min , this 
integral can be extended to the range — oo < x < oo. The operator D = d/dx is then 
seen to be skew Hermitian and «S + «S_ is Hermitian. 

Now, if the function ty^ S2M is sufficiently smooth, is non-zero over a sufficiently 
narrow region of x within the limits m m j„ < x < m max and is centered about a value 
Xq, we can make the so-called SHA [H [2| which, in addition to the continuous variable 
approximation, consists of dropping all terms in F, fo and f\ for which the expansion 
of the operator S + S- will be more than bilinear in x — Xq and d/dx. The conditions 
under which the SHA are valid are shown, in the following, to be well satisfied, for finite 
values of n, in the asymptotic limit as si — > oo and S2 — > 00. Thus, the SHA gives the 
asymptotic expression 

S + S_ l -A^- 2 + C(x - x ) - \b{x - x ) 2 + . . . , (40) 

2 dx z 2 

where 

E = F(x ), A = 2/i (x ), C = F'(x ), B = -F"(x ). (41) 

An examination of the values of these parameters reveals that A > and B > in 
situations of interest. Thus, we consider "H ~ — <S + <S^ and determine x to be the value 
for which C = 0, bringing "H to the form 

1 d 2 1 



n = -E + 



2a 2 dx 2 2 



-a 2 (x — x ) 2 



hu, (42) 



where 



(hu) 2 = AB, a 4 = j. (43) 

The eigenfunctions for this Hamiltonian are harmonic oscillator wave functions, and the 
eigenvalues are 

E n = -E + (n + \)hw. (44) 

Thus, we obtain the desired asymptotic Clebsch-Gordan coefficients in the form given by 
Eqn. (T5]). Note that the eigenfunctions of a Hamiltonian are only ever defined to within 
an arbitrary phase. There are also arbitrary phases in the definition of CG coefficients. 
The phases in Eqn. ([2l) were chosen to reproduce the standard phases of the Condon 
and Shortley CG coefficients [T9] . 



2.2. Simplified analytical expressions for asymptotic SU(2)CG coefficients 

The SHA results given above are easily calculated, and give accurate results in the 
limit as s\ — > 00 and s 2 — > 00 but M and n remain finite. Simpler analytic asymptotic 
expressions are obtained if we also neglect terms that go to zero as Si — > 00 and s 2 — > 00. 
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with 



In these limits 

Id 2 1 
H ->■ const. - 2 A °~^2 + C °( x ~ x °) + 2 B o( x ~ x °) 2 ' ( 45 ) 

M 2 \ 

A ~ 2a 1( T 2 1 - ■ , (46) 



(CTI + (T 2 ; 



C ~ 2 (^+^) M , (47) 

Bo „ 2 (* + ■»>' . (48) 
Thus, we can evaluate the asymptotic Clebsch-Gordan coefficients from Eqn. (j2J) with 

a 4 = ^0 = (gl + ^2) 4 , , 

i^O CTi + cr 2 
to obtain the expression 

(si,M—m, s 2 rra|s 1 + s 2 — n, M) 

~ ( " 1)n (7^!)^ n(a(m " Xo))e " |a2(m " xo)2 ' (51) 

and that of Eq. (JSJ) . 

Numerical results for su(2) CG coefficients 

We have ascertained and the following examples illustrate that the SU(2) CG coefficients 
satisfy the conditions for the validity of the SHA in the specified asymptotic limit as 
si, S2 — > oo for finite values of n and M. The following results show that in addition 
to being precise in these limits, the SHA and simplified SHA wave functions also give 
remarkably accurate values for modest values of s\ and s 2 and surprising large values of 
M. 

The following figures compare the values of exactly computed SU(2) CG coefficients, 
(si,M — m, s 2 m|si + s 2 — n,M), with the SHA expressions given by Eqn. (J2J). The 
continuous (red) lines are those of the full SHA approximation with the parameters as 
defined in Sect. 12.11 The dashed (black) lines are those of the simplified SHA coefficients 
with the parameters given by Eqns. (1491 - 1501) . 

Figure [1] illustrates that, for s\ and s 2 as small as 20 and 15, respectively, both 
the full and simplified SHA yield approximate SU(2) CG coefficients for n < 5 that are 
almost indistinguishable from each other. The figure shows that inaccuracies become 
visible for n = 5 coefficients when \m\ > 8. The results become increasingly accurate 
for larger values of s\ and s 2 and are precise in the asymptotic limit. 

Figure [2] illustrates that, for si = 60, s 2 = 40, and n = the simplified SHA is very 
accurate for \M\ < 30 whereas the full SHA is accurate for all the values of M shown. 
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n=0 




Figure 1. The Clcbsch-Gordan coefficients (20, — m, 15,m|35 — n, 0) shown as a 
function of m for three values of n. Exact values are shown as dots, full SHA values 
as continuous (red) lines, and simplified SHA values as dashed lines. 




m+41 

Figure 2. The Clebsch-Gordan coefficients (60, M — m, 40, m|100, M) shown as a 
function of m for a range of M values. Exact values are shown as dots, full SHA values 
as continuous (red) lines, and simplified SHA values as dashed lines. 

Figure |3] illustrates how accurate the SHA and simplified SHA Clebsch-Gordan 
coefficients can be for quite small values of s\ and S2 provided n and \M\ are kept even 
smaller. The region in which the results are at their worst is for values of m close to 
its boundary values, especially in situations in which the asymptotic expressions extend 
beyond these boundaries. 

3. Asymptotic SU(1,1) Clebsch-Gordan coefficients 

We now consider the operators {K+, K-, Ko}, that satisfy the commutation relations 
of the complex extension of the su(l,l) Lie algebra: 

[k , K ± ] = ±K ± , [£_, K + ] = 2K . (52) 
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n=0 




Figure 3. The Clebsch-Gordan coefficients (10, —m, 7, m\17—n, 0) shown as a function 
of m. Exact values are shown as full dots, full SHA values as continuous (red) lines, 
and simplified SHA values as dashed lines. 



Basis states, {\kn),n = 0, 1,2, . . .}, for a unitary su(l,l) irrep of the positive discrete 
series are defined by the equations 

K \kn) = (k + n)\kn), 

K + \kn) = ^{2k + n)(n + l) \k,n+l), 

K_\kn) = ^{2k + n-l)n \k, n - 1), (53) 
K + K_\kn) = (2k + n-l)n\kn). 

Basis states for irreps of su(l,l) within the tensor product of two such irreps are 
given by linear combinations 

\Kn) = 2^ 1^2) <8> \kini) C(k u ni, k 2 , n 2 ; K, n), (54) 

where C(ki,ni,k2,ri2;K,n) is an SU(1,1) CG coefficient in the notation of Van der 
Jeugt [13], equal to the overlap of a coupled and uncoupled tensor product states 

C(kx,n 1 ,k2,n 2 ; K,n) = [(A;i,ni| <g> (k 2 , n 2 \]\Kn). (55) 

Thus, if we put n 1 = |iV — m and n 2 = |iV + m and denote an uncoupled tensor product 
state by 

Ik^Nm) = \k 2 , IN + m) ® \N - m), (56) 

we obtain the more useful expression of a CG coefficient as an overlap 

C(ki,ni, k 2 ,n 2 ; K,n) = (kik 2 Nm\Kn) . (57) 

These coefficients are defined (to within arbitrary phase factors) by the requirement 
that the states {\Kn)} satisfy the eigenvalue equations 

K \Kn} = (K + n)\Kn), (58) 
K + k_\Kn) = (2K + n- l)n\Kn), (59) 
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with 



Ao = ^o + ^o. K ± = kl + k 2 ± 



(60) 



and the understanding that 



Kl(\k 2 m 2 ) <8> \k\mi)) = \k 2 m 2 ) <8> Kl\kiirii), 
Kl(\k 2 m 2 ) ® \kimf)) = K*\k 2 m 2 ) <g> \hrnf). 



(61) 
(62) 



3.1. SU(1,1) CG coefficients in the shifted harmonic approximation 

We now determine asymptotic limits to these CG coefficients in the SHA and show them 
to be precise for large values of N and finite values of n and \k\ — k 2 \. 

Before embarking on an SHA calculation of asymptotic SU(1,1) CG coefficients, 
we first examine some coefficients to see if they satisfy the necessary conditions for the 
validity of the SHA. From exact calculations in the phase convention of Van der Jeugt 
[20] , it is determined that, for large values of N and small values of n, the SU(1,1) CG 
coefficients (k\k 2 Nm\Kn) , when multiplied by the phase factor (— l) Ar+m and regarded 
as functions of m, approach standard harmonic oscillator wave functions. Thus, we 
define an overlap function of the discrete variable m by 



with K = N -n + ki + k 2 . 

Because the state \Kn) is an eigenstate of the operator K+K_, it follows that the 
representation of this state by the function i/j^ lk2N , of the discrete variable m is an 



,kxk 2 N 



(63) 



■i i 



eigenfunction of the operator /C + /C_ defined by 



/C + /C_^ lfc27V (m) = (-l) N ^ m {k 1 k 2 Nm\k + k^\Kn) 

= Y,(- 1 ) m ~ P ( k i k 2Nm\k + k.\k 1 k 2 Np)^ lk2N (p). (64) 



v 




but now with 



(kxkiNmlK+K^kxkiNm) 
(2k x + \N -m-l)(\N -m) 



fi{m) 



+ (2k 2 + \N + m - 1){\N + m), 

- (k l k 2 Nm\k + k_\kik 2 N,m + 1} 

- [{2k x + \N-m- 1){\N - m) 



(66) 



x (2k 2 + \N + m){\N + m + 1)] ' 



(67) 
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To determine asymptotic expressions for the functions, ip klk2N (m) , as eigenf unctions of 
/C + /C_, we extend these functions of the discrete variable m to functions, vJ/^i^a^ of a 
continuous variable x with the property that 

m k n ^ N (x) = ^ k n lk2N {x), (68) 

whenever x is in the domain of the discrete variable m. Thus, as for SU(2), we obtain 
an expression for /C + /C_ as the differential operator 

/C+/CL = F(x) + Dfi(x)D, (69) 
where F(x) = fo(x) + fi(x) + fi(x — 1) and D = d/dx. With k± and k 2 defined by 

2«i := 2k x + |iV, 2k 2 := 2k 2 + ~iV, (70) 
we also obtain the expressions 

f (x) = (2«i - x - l)(|iV - x) + (2k 2 + x - l)(|iV + x) (71) 

h{x) = - [{2 Kl -x-l){\N- x){2k 2 + x){\N + x + 1)]'. (72) 

Provided the extension of ip klh2N (m) to the smooth function \E^ lfc2Ar (x) does not 
require the latter to be non-zero for any x that is outside of the limits for m, it is seen 
that fi(x) is real for all x for which \E^ lS2M (x) is non-zero. The limits on the values 
of m are seen, from Eqn. ( 1561) . to be such that —N/2 < m < N/2. Then, because the 
norm of the function ip klk2N (m) is given by 

N/2 

ll^ lfc27 V)|| 2 = E l^ lfc27 V)| 2 , (73) 
m=-N/2 

it follows that the corresponding smooth function ^/ h l lk2N (x) should have norm given by 

r N/2 

||tf*ik*||2 = / ^^(x)] 2 ^. (74) 



J —N/2 

Also, when evaluated without approximation, ^/ k l lk2N (x) is zero for all x > N/2 and all 
x < —N/2 and this integral can be extended to the range — oo < x < oo. The operator 
D = d/dx is then seen to be skew Hermitian and /C + /C_ is Hermitian. 

Now, if the function ty k l lk2N (x) is sufficiently smooth, is non-zero over a narrow 
region of x within the limits —N/2 < x < N/2, and is centered about a value xq, we 
can again invoke the SHA of dropping all terms that are more than bilinear in x — Xq 
and d/dx in an expansion of the operator /C + /C_. This gives 

/C + /C_ ^E- \a^ 2 + C(x - x ) + l -B(x - x ) 2 , (75) 

where 

E = F(x ), A — — 2/i(x ), C = F'(x ), B = F"(x ). (76) 

This approximation becomes precise, in the N — > oo asymptotic limit, provided the 
parameters n, k\, and h 2 remain finite. In fact, it turns out that h\ + k 2 can also be 
large provided the difference \k\ — k 2 \ remains small in comparison to N. The manner 
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in which the SHA ceases to be valid, for large values of n and \k\ — k 2 \, relative to N, 
is shown in Sect. 13.31 Thus, we consider the SHA Hamiltonian 

H = E- \a^- 2 + C(x - s ) + l -B(x - x ) 2 . (77) 

The appropriate value for x® is that for which C = and 

1 d 2 1 



U = E 
with 



2a 2 dx 2 2 



+ -a 2 (x — xq) 2 



fiw, (78) 



(hoo) 2 = AB, a 4 = j. (79) 

The eigenfunctions for this Hamiltonian are harmonic oscillator wave functions, and the 
eigenvalues are given by 

E n = E + (n + \)froj. (80) 

With a standard choice of phase for the harmonic oscillator eigenfunctions, we then 
obtain the desired asymptotic Clebsch-Gordan coefficients, as given in Eqn. flSJ). 

3.2. Simplified analytical expressions for asymptotic SU( 1,1) CG coefficients 

The SHA results given above are easily calculated, and give accurate results for large 
values of N — > oo and small values of n. Simpler analytical expressions are obtained 
if we further neglect terms in the SHA Hamiltonian which go to zero in the N — > oo 
asymptotic limit. In these limits 

1 d 2 1 

U ->■ const. A^-— + C x + -Bx 2 , (81) 

2 dx 1 2 



with 



A = 2N V ^, (82) 
+ 



C = — Z- — ( K 2 ~ Ki). (83) 



Thus, we can evaluate the asymptotic Clebsch-Gordan coefficients from Eqn. (jHJ) with 



a 2 



\J Bq/Aq and xq = —Co/B given explicitly in Eqn. 



3.3. Numerical results for su( 1,1) CG coefficients 

We have ascertained and the following examples illustrate that the SU(1,1) CG 
coefficients, when multiplied by a phase factor ( — l) N+m , satisfy the conditions for the 
validity of the SHA in the specified asymptotic limit as N — > oo for finite values of 
n and \k\ — A:2 1 - The restriction on the value of n is because the harmonic oscillator 
wave functions become broad and oscillate rapidly for large values of n with the result 
that if n/N is too large the conditions for the validity of the SHA cease to be satisfied. 



SHA and asymptotic SU(2) and SU(1,1) CG coefficients 



14 



Likewise, the restriction on the value of \k\ — k 2 \ is understood to arise from the value 
of Cq which, for iV — > 00, approaches C ~ 3(k 2 — kf). Thus, unless \ki — k 2 \ « N, 
it is not guaranteed that the centroid, xq ~ —Cq/Bq, of the SHA wave function will be 
sufficiently close to the center of the domain —N/2 < x < N/2 for the SHA to be a valid 
approximation. However, the results below show that in addition to being precise in 
the asymptotic limits, the SHA (as opposed to the simplified SHA) remains remarkably 
accurate for relatively large values of \k\ — k 2 \ when N is large. 

The following figures compare the values of exactly computed SU(1,1) CG 
coefficients, (kik 2 Nm\Kn), with the SHA expressions given by Eqn. ([8]). The continuous 
(red) lines are those of the full SHA approximation with the parameters as defined in 
sect. 13 . 11 The (black) dashed lines are those of the simplified analytical SHA coefficients 
with parameters given by Eqn. (Q. The asymptotic limits require that n remains finite 
while N — y 00. The simplified SHA coefficients retain the property of being precise in 
the N — y 00 limit but requires, in addition, that \k\ — k 2 \ remains small compared to 
N/2. In fact, as the following figures illustrate the SHA coefficients yield surprisingly 
good approximations for quite modest values of N and the analytical approximations 
are seen to be almost indistinguishable from the full SHA expressions when \k\ — k 2 \ is 
not too large. 

Figure |4] illustrates the accuracy that can be obtained with the above-defined 
asymptotic SU(1,1) CG coefficients for N = 100, n — 10 and k 2 — k\ = 7. It 
shows that the SHA coefficients and the analytical approximations to them are virtually 
indistinguishable. It also shows that the errors in the asymptotic coefficients for a large 
but finite value of N start to become most noticeable, for these relatively large values 
of n and \ki — k 2 \, at the upper and lower reaches of m. 



-0.1 



0.2 



0.1 







10 




90 
m+51 



100 



Figure 4. SU(1,1) Clebsch-Gordan coefficients (-l) m (10, 17, 100, m|117, 10) shown as 
a function of m. Exact values are shown as full dots, full SHA values as continuous 
(red) lines, and simplified SHA values as dashed lines. 



Figure shows the SU(1,1) CG coefficients for a range of ki — k 2 values. The 
calculations for other k\ — k 2 values show that the simplified SHA coefficients are close 
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to those of the full SHA for \k\ — k%\ < 25 but are noticeably different for larger values 
of \k\ — k 2 \ as seen for the (fci, k 2 ) = (5, 50) and (50,5) coefficients shown in the figure. 
It is also seen that, for the full SHA, errors start to become evident for large values of 
\k\ — k 2 \ as the value of m approaches its upper or lower bounds. 




10 20 30 40 50 60 70 80 

m+41 



Figure 5. SU(1,1) Clebsch-Gordan coefficients (-l) m (fcifc 2 , 80, m\K, 0) with K = 
80 + ki + &2 for a range of (^1,^2) values shown as functions of m. Exact values 
are shown as full dots, full SHA values as continuous (red) lines, and simplified SHA 
values as dashed lines. The simplified SHA values for \ki — > 45, are unacceptably 
inaccurate and are not shown. 

The SU(1,1) CG coefficients, (kik2Nm\Kn) , are given exactly both in the SHA 
and in the simplified analytical approximation to the SHA, for finite values of n and 
\k\ — k 2 \, in the iV — > oo asymptotic limit. However, the full SHA expressions remain 
noticeably more accurate over a considerably larger domain. In fact, the coefficients 
continue to be accurate for relatively small values of N and even smaller values of n and 
\k\ — k 2 \- This is illustrated for N = 10 in fig. [6j The region in which the results are at 



n=0 




Figure 6. SU(1,1) Clebsch-Gordan coefficients (-l) m (l/2, 3/2, 10, m|12, 0) and 
(1/2, 3/2, 10, m|10, 2) shown as a function of m. Exact values are shown as full dots, 
full SHA values as continuous (red) lines, and analytical asymptotic limits of the SHA 
values as dashed lines. 



their worst is for values of m close to their limits, especially in situations in which the 
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asymptotic expressions extend beyond these limits. 



4. Comparisons with the random-phase approximation (RPA) 

Whereas the SHA and RPA are both harmonic oscillator approximations and both 
become precise in asymptotic limits, they apply in different complementary situations. 
In their application to the derivation of asymptotic SU(2) and SU(1,1) CG coeficients, 
the RPA can be viewed as a contraction of the su(2) Lie algebra to a harmonic 
oscillator boson algebra. Thus, we consider the possibility that the SHA might also 
correspond to such a contraction, albeit one that is valid in a different domain of an 
SU(2) representation space. 

4.1. SU(2) CG coefficients 

We start from the observation that SU(2) CG coefficients are the eigenstates of a 
Hamiltonian 

H = aS + xS+S^ (85) 

on the tensor product space of two su(2) irreps, {si} ® {$2} in a basis of product states 
{\sirrii) <g) \ s2rri2)}. The Hamiltonian H has eigenstates, {\SM)}, with eigenvalues given 
by 

E SM = aM + X [S(S + l)-M{M-l)]. (86) 

Inspection shows that, when x > ~~ a i f° r a > 0, the lowest-energy eigenstate of H is 
the state with S — a, M = —S, where a = s\ + S2 and, when x < ~ct, the lowest-energy 
eigenstate is the state with S = a, M = 0. Thus, the two situations are relevant for the 
calculation of asymptotic SU(2) CG coefficients {s\m\ S2m2\SM), for large values of S\, 
S2, and S, with M close to either —S or 0, respectively. As we now observe, RPA gives 
solutions in the first scenario whereas the SHA does so in the second. 

In the a — > 00 asymptotic limit for a > and x > ~ a , h is convenient to relabel 
the low-lying states by \SM) — > \nm), where n = a — S and m = S + M. The energy 
£ nm = Esm is then given to leading order in n and m by 

£nm ~ —ac + an + [a + x(2o" + (87) 

which is an eigenvalue of a Hamiltonian expressed in terms of the raising and lowering 
operators of two simple harmonic oscillators by 

U = -aa + ac ! c+ [a + x(2<r+ l)]6 f 6, (88) 

with 

cc\nm) = n|nm), b^b\nm) = m\nm). (89) 

For asymptotically large values of Si and s 2 , this result can be obtained in the RPA by 
contracting the su(2) algebra in each of the {s\} and {^2} irreps by 

S + ~ a\ SL ~ —= a, S = —s + a) a, (90) 
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for s = Si and s = s 2 , respectively. Diagonalization of the Hamiltonian H of Eqn. (1851) 
in this contraction limit, then leads to the asymptotic expression f|88|) and determines 
the c and b boson operators from which one can derive asymptotic CG coefficients. The 
contraction (190]) which leads to this RPA result is valid in the domain of states, {\SM)}, 
in which M is close to —S and S is close to si + s 2 , for large values of si and s 2 . 

For x < ~ a i the RPA breaks down because the ground state of the Hamiltonian 
suddenly flips from an \S,M — —S) state to an \S,M = 0) state at x — ~ a - However, 
for x < ~~ a the SHA provides asymptotic solutions. 

In the a — > 00 asymptotic limit for x < —ct-, if is convenient to relabel the low-lying 
states by \SM) — > \nM), where n = o — S . It is then found that S n M = E$m is given 
to leading order in n by 

Sum w X°(° + 1) - X(2a + l)n + aM - X M(M - 1). (91) 

Thus, in this asymptotic limit, the spectrum of eigenvalues of the Hamiltonian (185|) is 
that of a simple harmonic oscillator coupled to a U(l) rotor 

H = x°(° + 1) - X(2a + l)Jc + a^o - xSo(S - 1), (92) 

where 

c ] c\nM) = n\nM), S \nM) = M\nM). (93) 

The SHA gives an explicit expression for a shifted harmonic oscillator that is 
essentially equivalent to %. We have not succeeded in deriving this Hamiltonian by 
a contraction of the su(2) algebra. However, it is noted that the SHA is based on a 
realisation of the su(2) Lie algebra given by the expressions, in terms of the harmonic 



oscillator operators, x = x, p = —ihd/dx 

S + ^S+ = e~ i?5 [s(s + 1) - x{x + 1)] 1/2 , (94) 

5_ -> 5_ = [s(s + 1) - x{x + l)] 1/2 e ip , (95) 

Sq^Sq = x. (96) 



Thus, it would appear likely that the SHA Hamiltonian (!9~2l could be obtained as a 
contraction of this realisation. 

4.2. SU(1,1)CG coefficients 

We now compare the complementary RPA and SHA derivations of asymptotic SU(1,1) 
CG coefficients by diagonalization of a Hamiltonian 

H = ak + X K+K-. (97) 

The spectrum of eigenstates of this Hamiltonian are the SU(1,1) states {\Kn)} with 
eigenvalues given by 

E Kn = a(K + n) + X (2K + n - l)n, (98) 

or, with a relabelling of states by \Kn) — > \Nn), where N = K + n — ki — k 2 , by 

£nu = a(h + k 2 + N) + X (2h + 2k 2 + 2N-n- l)n. (99) 
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For a and x positive, the lowest value of S^ n = E^ n is for n = N = 0. Then, for 
asymptotically large values of k\ and k 2 , the low-lying eigenvalues are given to leading 
order in N and n by 

£ Nn ~ a(h + k 2 ) + aN + X (2h + 2k 2 - l)n, (100) 
which are the eigenvalues of the Hamiltonian 

H ~ a(h + k 2 ) + ac f c + x(2Jfei + 2£; 2 - l)¥b (101) 
for a two harmonic oscillator system with 

Jc\Nn) = N\Nn), tfb\Nn) = n\Nn). (102) 

The corresponding eigenvectors, in the space of coupled tensor product states 
{\kinf) (g) \k 2 n 2 )}, and hence a subset of asymptotic SU(1,1) CG coefficients, can be 
obtained in the RPA by a contraction of the two su(l,l) irreps {k%} and {k 2 }. The 
relevant contraction is obtained by observing that for small values of n, the su(l,l) 
relationships 

[K_,K + ]\kn) = 2K \kn) = 2(k + n)\kn), [S ,§±] = ±S±, (103) 
are satisfied in the k — > oo asymptotic limit, for finite values of n, by 

K + ~ a\ —= K_ ~ a, K ~k + a f a. (104) 



'2k V2k 

Thus, by expressing the Hamiltonian of Eqn. fl97|) in terms of the su(l,l) operators 
of the {ki} and {k 2 } sub-representations, it becomes the Hamiltonian for two coupled 
harmonic oscillators which, in the RPA, are decoupled by diagonalization to give the 
uncoupled oscillator expression of Eqn. (llOip . The contraction fll04p which leads to this 
RPA result is valid in the domain of coupled states, {\Kn)}, in which K is close to 
k\ + k 2 and n is small. 

Other su(l,l) CG coefficients are obtained by considering the eigenvectors of the 
Hamiltonian (I97p in the subspace of states of a fixed value of N. For asymptotically 
large N and small n, Eqn. ( 199|) then reduces to the eigenvalues 

£nu ~ a(h + k 2 + N)+ 2 X (h + k 2 + N)n (105) 

of the Hamiltonian 

H ~ a(h + k 2 + N) + 2 X (h + k 2 + N)b ] b (106) 

for a simple harmonic oscillator with b'b\Nn) = n\Nn). The SHA gives an explicit 
expression for this Hamiltonian as a shifted harmonic oscillator Hamiltonian. We have 
not succeeded in deriving this expression in terms of a contraction of the su(l,l) algebra. 
However, the SHA is based on an approximation to a realisation of the su(l,l) Lie 
algebra given by the expressions, in terms of the harmonic oscillator operators x = x, 
p = —ifrd/dx, 

K + ^t+ = e~ ip [{2k + x){x + 1)] 1/2 , (107) 
K_ -)• t- = {(2k + x){x + l)] 1/2 e ip , (108) 
K ^]C =k + x. (109) 
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Thus, it would appear once again likely that the Hamiltonian f )106p can be obtained as 
a contraction of this realisation. 

5. Discussion and conclusion 

In this paper we have demonstrated that turning a three-term recursion relation into 
a second order differential equation makes it easy to understand the oscillatory nature 
of CG coefficients. From a more numerical perspective, it has been noted that the 
final forms given in the paper, even the simplified expressions, are remarkably accurate 
even for values of the parameters that are far from the asymptotic limits in which they 
become precise. 

It is readily ascertained that our asymptotic SU(2) CG coefficients retains the 
symmetry 

(si^i s 2 m 2 \s 1 -\-s 2 - n, M) = (-l) n (s 2 m 2 Sim 2 |si + s 2 - n, M) (110) 

under the exchange s±mi -H- s 2 m 2 . Indeed, this is seen from the simplified expression of 
Eq. (j3J) in which the argument of the Hermite polynomial simply changes its sign under 
this exchange. A parallel result is obtained for an SU(1,1) CG coefficient for which, in 
the sign convention used, 

C(ki,n 1 ,k 2 ,n 2 ; K,n) = {-l) n C{k 2 ,n 2 ,k 1 ,n 1 ]K,n). (Ill) 

Due to the asymmetrical way in which the asymptotic limits are approached, it is not 
expected that the more general symmetries of these CG coefficients will be preserved. 
However, such symmetries as are known, e.g., for the exchange s\mi <H- S3, — 7713, can be 
used to re-arrange the arguments of a given CG coefficient such that the value of n is 
minimised. The effective range over which the asymptotic approximations are expected 
to produce acceptable results is thereby increased. 

More generally, the successes of the SHA in its applications to date [H [2j [3] 
suggest that it is a potentially powerful technique that could be applied more generally. 
Note, for example, that the finite and discrete series of irreps of all semi-simple Lie 
algebras are characterized by the irreps of their many su(2) and/or su(l,l) subalgebras 
generated by raising and lowering operators. In particular, every pair of raising and 
lowering operators, X^ 1 , of a semi-simple Lie algebra can be normalized to have SU(2) 
commutation relations [X+,X~] = h v , [h^jX^] = ±2X^. Note also that many 
dynamical systems have low-energy collective states that are well approximated as 
harmonic vibrational states. As a result the RPA (random phase approximation) 
has become a powerful tool in many-body theory. It is also of interest to explore 
the possibility of mapping the algebraic structure of a problem to a rotor algebra, for 
example, rather than that of a harmonic vibrator. 

The search for an extension of the SHA to apply to many-body systems with more 
general algebraic structures is worthwhile because it is known that, for a variety of many- 
body systems, the RPA works well for describing vibrational normal mode excitations 
with relatively weak interactions but breaks down when the interactions become too 
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strongly attractive. This is the situation in which the system is understood to undergo 
a phase change; the frequency of one of the vibrational modes goes to zero and a 
deformed equilibrium state, with rotational plus vibrational degrees of freedom, emerges. 
Thus, we are optimistic that a generalised SHA will prove to be appropriate for such 
situations. Section 4 has shown the SHA to be complementary to the RPA its ability to 
give (asymptotic) solutions to an eigenvalue equation in regions where the RPA breaks 
down. Moreover, the SHA has been shown, in previous applications, to provide practical 
solutions of the multi-level BCS Hamiltonian that are remarkably accurate for relatively 
strong pairing interactions [3]. 



1 f s\ + s 2 



Appendix A. A special case 

The special case 

(si, -to; s 2 , mlsi + s 2 , 0) 

\7T \ SiS 2 

follows by careful application of Stirling's formula 

to the exact expression 

(si, -to, s 2 ,m\si + s 2 ,0) 

(2 Sl )!(2s 2 )!( Sl + S2 )!( Sl + S2 )! 



-m 2 (si+S2)/2s\S2 



SI 




s s+1 e~ s 



_(2si + 2s 2 )!(si + m)!(si - m)!(s 2 + rn)\(s 2 - to)! 
First, one can use Eqn. (IA.2l) to show that 

„2s+2 



S\S\ 



[s + my.is — my. 



s 2 — m? 



[s + m) s+m+1 (s - m) s - m+1 



With x = m/s, this expression reduces to 



s\s\ 



(s + m)\(s — m)\ y/Y 



x' 



(l-x) x 



1 + xY 1 -x 2 



and, for a finite value of m, reduces further in the limit as s — > oo to 



s\s\ 



-m 2 /s 



(s + m)\(s — m)\ (l+m 2 /s 2 ) s 
Eqn.( 1A.3l) can then be manipulated to directly yield Eqn. (1A.ljl 
As a special case we observe that 

(2s)\(2s)\ 



[s, —m; s, m\2s, 0) 



a/(4s)! (s + m)\(s - to)! 



ITS 



-m 2 /s 



Both Eqn.f lA.l]) and Eqn. flA.8j) agree with Eqn.(T2]) for n = and M 
limit where s±, s 2 — > oo. 



(A.l) 



(A.2) 



. (A.3) 



(A.4) 



(A.5) 



(A.6) 



(A.7) 

(A.8) 
in the 
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